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1. Introduction 



It is an essential step for lattice QCD calculation to reproduce the harden spectrum with vari- 
ous systematic errors under full control and hence show that QCD is the fundamental theory of the 
strong interaction. Recent developments in simulation algorithms and computational facilities now 
allow us to remove the most troublesome systematic errors: the quenching effect and chiral extrap- 
olation. The main task of the previous CP-PACS/JLQCD project [[I], ^ was to remove the former by 
performing A/y = 2 + 1 lattice QCD simulations with the nonperturbatively 0(fl;)-improved Wilson 
quark action and the Iwasaki gauge action on a (2fm)^ lattice at three lattice spacings. The 
lightest up-down quark mass, however, was still 64MeV corresponding to nij^/mp ^ 0.6 so that 
there remained an ambiguity from a long chiral extrapolation. 

In the PACS-CS (Parallel Array Computer System for Computational Science) project |^ ^, 0, 
|8]], we aim to perform Nf = 2 + 1 lattice QCD simulations at the physical point using the PACS-CS 
computer with a total peak speed of 14.3 TFLOPS developed and installed at University of Tsukuba 
on 1 July 2006. The quark and gauge actions are the same as in the previous CP-PACS/JLQCD 
project. To reduce the degenerate up and down quark masses we employ the domain-decomposed 
HMC (DDHMC) algorithm with the replay trick proposed by Liischer |^ |l^], whose effectiveness 
for small quark mass region has already been shown in the Nf = 2 case |lT]] incorporating the 
multiple time scale integration scheme [p^]. The strange quark is included by the UV-filtered 
Polynomial Hybrid Monte Carlo (UV-PHMC) algorithm [Q. 

In this report we present simulation details and some preliminary results for the hadron spec- 
trum. The analysis on the chiral behavior of the pseudoscalar meson masses and the decay constants 
using chiral perturbation theory is given in a separate report [14]. 



2. Simulation details 



We employ the 0(a) -improved Wilson quark action with a nonperturbative improvement co- 
efficient csw = l-715[|l5l] and the Iwasaki gauge action at j8 = 1.90 on a 32^ x 64 lattice. Table |] 
summarizes our simulation parameters. Simulations are carried out for six combinations of the 
hopping parameters (K"ud,K's). The heaviest case of (KTudjKTs) = (0.13700,0.13640) is chosen for 
a direct comparison of the PACS-CS results with the previous CP-PACS/JLQCD ones which em- 
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Table 1: Summary of simulation parameters. MD time is the number of trajectories multiplied by the 
trajectory length T. Tint[P] denotes the integrated autocorrelation time for the plaquette. 
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ployed this parameter as the smallest quark mass. Among two choices of K^, = 0.13640 is 
the physical point = 0.136412(50) as determined in the previous CP-PACS/JLQCD workl]!), H], 
while K"s = 0.13660 is to investigate the strange quark mass dependence. 

The DDHMC algorithm is implemented by domain-decomposing the full lattice with a 8^ 
block size as a preconditioner for HMC. The domain-decomposition factorizes the up-down quark 
determinant to the UV and the IR parts distinctively. On the other hand, the UV-PHMC algo- 
rithm for the strange quark is not domain-decomposed. We can incorporate the multiple time scale 
integration scheme in these algorithms to reduce the simulation cost efficiently. The relative mag- 
nitudes of the force terms are found to be 

||iv||:||Fuv||:||fiR||~ 16:4:1, (2.1) 

where Fg denotes the gauge part and Fuv,ir for the UV and the IR parts of the up-down quarks. We 
write the associated step sizes for the above forces as 5Tg = t/NqNiNj, 5tuv = 5xir = 

Z/N2 with z the trajectory length, where the integers No,N\ ,N2 are chosen such that 

5Tg||Fg|| ^ 5tuv||Fuv|| « 5tir[|Fir||. (2.2) 

We take No = N\ = 4. Since we have found ||Fs[| ||fiR|| for the strange quark part, we choose 
5Ts = 5tir. The value of N2 and the polynomial order for UV-PHMC A'poiy are adjusted taking 
account of acceptance rate and simulation stability. 

The inversion of the Wilson-Dirac operator D is carried out by the GCR solver with the stop- 
ping condition |D;c — < 10^^ for the force calculation and 10^^^ for the Hamiltonian, which 
guarantees the reversibility of the molecular dynamics trajectories to high precision: \AU\ < 10^^^ 
for the link variables and \AH\ < 10^^ for the Hamiltonian. 

3. Plaquette history and autocorrelation time 

In Fig. [1] we show the plaquette history and the normalized autocorrelation function at {K^d, K"s) = 
(0.13727,0.13640) as a representative case. The integrated autocorrelation time is estimated as 




500 1000 1500 2000 2500 "°-^0 50 100 150 



Figure 1: Plaquette history (left) and normalized autocorrelation function (right) for (K'ud;K's) = 
(0.13727,0.13640). 
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Figure 2: {am,^)^ (left) and amp (right) as a function of 1 / K"u(j. Black and red symbols denote the PACS-CS 
and the CP-PACS/JLQCD results, respectively. 





lattice size 


aniji 


amp 




PACS-CS 


32^ X 64 


0.3220(6) 


0.506(2) 


0.726(3) 


CP-PACS/JLQCD 


20^ X 40 


0.3218(8) 


0.516(3) 


0.733(4) 



Table 2: PACS-CS and CP-PACS/JLQCD results for it, p and nucleon masses at (K'ud,K's) = 
(0.13700,0.13640). 



'^mAP] = 20.9(10.2) following the definition in Ref. For other hopping parameters we have 
found similar behaviors for the normalized autocorrelation function. Although we hardly observe 
the quark mass dependence for Tint [P] listed in Table the statistics may not be sufficiently large 
to derive a definite conclusion. 



4. Hadron spectrum 

We measure hadron correlators at every 10 trajectories at the unitary points where the valence 
quark masses are equal to the sea quark masses. Light hadron masses are extracted from single 
exponential fits to the correlators with an exponentially smeared source and a local sink. Sta- 
tistical errors are estimated by the jackknife method, whose bin size is chosen to be 50 molecular 
dynamics time based on the bin size dependence of the errors. 

We first compare the PACS-CS results on 32^ x 64 with the previous CP-PACS/JLQCD results 
on 20^ X 40||l], ^ at j8 = 1.90. Figure ^ shows the up-down quark mass dependence of (aniji)^ and 
amp with fixed at 0.13640. As for the pion mass we observe that the PACS-CS and the CP- 
PACS/JLQCD results are smoothly connected as a function of 1 / KTud- Table ^ shows that the results 
for the pion mass at KTud = 0.13700 are consistent within the error bars. On the other hand, we find 
1—2% deviation for the p meson and nucleon masses, which could be finite size effects. 

Figure ^ shows the hadron effective masses at K"ud = 0.13727 and 0.13781. We observe clear 
plateau for the mesons except for the p meson and also good signal for the baryons thanks to a large 
volume. Especially, the Q. baryon has a stable signal and a weak up-down quark mass dependence 
for our simulation parameters. Taking advantage of this virtue we employ the Q. baryon as input to 
determine the lattice cutoff. Combined with the additional inputs of mji and m^ to determine the 
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t t 

Figure 3: Effective masses for the mesons (top) and the baryons (bottom) at K^^ = 0.13727 (left) and 
0.13781 (right). 
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Figure 4: Light hadron spectrum extrapolated at the physical point with il-input (red) and ^-input (blue). 

physical up-down and strange quark masses, we obtain a^^ = 2.256(80)GeV. With the use of this 
cutoff we find that the lightest pseudoscalar meson mass we have reached is about 210MeV. The 
light hadron spectrum extrapolated to the physical point is shown in Fig. ^ and compared with the 
experimental values. They show a good agreement within the error bars, though our results contain 
possible 0((fl'AQCD)^) cutoff errors. In Fig. ^ we also plot the results for the 0-input case, which 
are consistent with the H-input case. A more detailed description about the determination of a^^ 
and the chiral extrapolation of the hadron masses is given in Ref. [p^]. 

We calculate the bare quark masses using the axial vector Ward-Takahashi identity (AWI) 
defined by am^*^ = hm,^^ (V4A^"'P(f)/'(0)) /{2{P{t)P{0))), where A™p is the nonperturbatively 
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Figure 5: Energy levels of the p meson and the two pion states with the momentum 27r/L as a function of 
the up-down quark mass with Ks = 0. 13640 (left) and time dependence of the R function (right). 



0(a) -improved axial vector current [16|. Employing the perturbative renormalization factors Za^p 



evaluated up to one-loop level [17, 18], we obtain 



m 



MS 
ud 



(At =2GeV) =2.3(1 l)MeV, ^(^t = 2GeV) = 69.1(25)MeV 



(4.1) 



The physical up-down quark mass is about half of our lightest one m^^{}X = 2GeV) = 5.61(40) 
at (jCud, KTs) = (0.13781,0.13640). We also measure the pseudoscalar meson decay constants, for 
which we find 



= 144(6)MeV, /^^ = 175(6)MeV, /^/^ = 1.219(22) (4.2) 

at the physical point with the perturbative Z^. They are about 10% larger than the experimental 
values, though their ratio, which is free from the ambiguities of the renormalization factors, is con- 
sistent with the experimental value. For both the quark masses and the pseudoscalar meson decay 
constants, our concern is the use of the perturbative renormalization factors which might cause 
sizable systematic errors. We are now calculating the nonperturbative Za^p with the Schrodinger 
functional scheme. 



5. p-TiTi mixing 

Since our simulations are carried out at sufficiently small up-down quark masses, it would 
be interesting to investigate the p-nn mixing effects [ [l9| , ^ ]. We find that, while the rest mass 
nip is always smaller than the two-pion energy 2\/m\ + {In/LY for all the hopping parameters, 
the energy of a moving p with a unit of momentum, i.e., \Jm^p + {2n/LY, becomes larger than 
\/m\ + (27r/L)2 + toward the lighter up-down quark masses. This situation is illustrated in 
Fig. I 

Let us consider two types of the p meson propagator with the momentum 2n/L: P||(27r/L) 
with polarization parallel to the spatial momentum and p±{2n/L) with polarization perpendicular 
to the spatial momentum. Phenomenologically the p-7i7i couphng is described by gpK7i£abcP^7l"d^n'^, 
which favors p|| (271 /L) 7i{2k/L)k{0) top±{2n/L) n{2n/L)n{0). We expect that the P||(27r/L) 
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propagator is more strongly affected by the mixing effects than the p±{27i/L) correlator. Since the 
mixing effects push up the upper energy level further and push down the lower energy level, they 
could be detected by measuring the R function defined by 

(Pii(/v)Pi;(p-o)> i.^.^^-,.,-.,,,, 

(p±(p,<)pI(p.o)) 

In Fig. ^ we plot log \R{t)\ as a function of t. The data show clear positive slopes which indicate 
fipil < Ep^^ . We also observe that the magnitude of the energy difference is rather small for JCud < 
0.13754, while it grows rapidly as the up-down quark mass is reduced for K^d > 0.13754. This 
feature may suggest that the (p|| (/7,f)P||(^,0)) correlator is getting dominated by the KK state 
toward the smaller up-down quark masses. In order to obtain a definite conclusion, we need more 
detailed investigations with increased statistics. 
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